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Asymptotic construction of pulses in the discrete Hodgkin-Huxley model for 
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A quantitative description of pulses and wave trains in the spatially discrete Hodgkin-Huxley 
model for myelinated nerves is given. Predictions of the shape and speed of the waves and the 
thresholds for propagation failure are obtained. Our asymptotic predictions agree quite well with 
numerical solutions of the model and describe wave patterns generated by repeated firing at a 
boundary. 

PACS numbers: 87.19.La,05.45.-a,82.40.Ck,02.30.Ks 



loa pne- 

m 



I. INTRODUCTION 

Understanding wave propagation in discrete excitable 
media is challenging because of poorly understood phe- 
nomena associated with spatial discreteness 1 , 2 ,P 
The study of the transmission of nerve impulses along 
myelinated axons is a paradigmatic example. Myelinated 
nerve fibers, such as the motor axons of vertebrates, are 
covered almost everywhere by a thick insulating coat of 
myelin. Only a fraction of the active membrane is ex- 
posed, at small active nodes called Ranvier nodes. The 
myelinated axons of motor nerves can be very long, and 
contain hundreds or thousands of nodes 6j. The wave 
activity jumps from one node to the next one giving rise 
to "saltatory" propagation of impulses J} . Saltatory con- 
duction on myelinated nerve models has two important 
features. One is the possibility of increasing the speed of 
the nerve impulse while decreasing the diameter of the 



nerve fiber 



. The other is propagation failure when the 



myelin coat is damaged 
multiple sclerosis. 



which causes diseases such as 



The propagation of nerve impulses along a myelinated 
fiber can be described by the spatially discrete Hodgkin- 

on 

Huxley system (HH) |3, llO|. For typical experimental 
data, this system couples equations for two fast vari- 
ables and two slow variables. Most analytical studies 
of action potentials have focused on discrete FitzHugh- 
Nagumo (FHN) models for one fast and one slow vari- 




ables 



for the leading edge 



and discrete bistable equations 



Careful computa- 



tional studies of models including more biological detail 
were carried out in 



lueis mciuu 

i], i3 3 



Numerical simulations show that reduced models in- 
volving only one fast and one slow variables are quanti- 
tatively inaccurate. Fig. ^ compares pulse solutions of 
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the full discrete Hodgkin-Huxley model (circles) and a 
FitzHugh-Nagumo type reduction (asterisks) generated 
by exciting the left end of a fiber. The pulse solutions 
of the HH model are slower and narrower than the pulse 
solutions of the FHN-like reduction. Discarding one of 
the two fast variables produces a pulse with an increased 
speed, as illustrated by Figure ^ To describe the prop- 
agation of the pulse leading edge we need to keep the 
two fast variables. Their evolution is described by a dis- 
crete bistable equation coupled to an ordinary differential 
equation. This system yields a Nagumo type equation (a 
single discrete bistable equation) only in a very particular 
limit. Similarly, discarding one of the slow variables pro- 
duces a slightly wider pulse than the one described by the 
two original slow variables. Moreover, FHN pulses and 
HH pulses have different structures. Pulse solutions of 
FHN-like models typically consist of two rigidly moving 
wave fronts ^ . HH pulses are "triangular waves" : they 
are formed by a leading wave front followed by a smooth 
region, as shown in Figure ^ 

In this paper we introduce an asymptotic strategy to 
construct solitary pulses and wave trains in the discrete 
HH model. Our asymptotic study exploits time scale sep- 
aration to split the variables in two blocks. The leading 
edge of the pulses is a wave front solution of the reduced 
system involving the two fast variables. This selects the 
speed of the wave. The two slow variables become rele- 
vant to determine the width of the peaks. Our asymp- 
totic constructions agree reasonably well with numeri- 




FIG. 1: Comparison of pulses for the HH model 
(circles), HH with + Uk = h* + n* (crosses), HH with 
(vk) (triangles) and the FHN-like system ©-0 
(asterisks): (a) D = 0.09, (b) D = 0.01. The pulses have been 
generated by a short stimulation at fc = and have been su- 
perimposed at a fixed time t. The pulse advances faster when 
only one fast variable is kept (asterisks and triangles). All 
variables are expressed in dimensionless units. 



cal solutions for typical experimental data. For contin- 
uous HH models, a quantitative approximation scheme 
exploiting time scale separation was proposed in Q. In 
the discrete case, traveling impulses do not appear as so- 
lutions of ordinary differential equations. Instead, more 
complicated differential-difference equations have to be 
analyzed. 

The paper is organized as follows. In Section II we de- 
scribe the discrete Hodgkin-Huxley model for myelinated 
nerves. We present a few numerical solutions and dis- 
cuss the reasons for the poor performance of FitzHugh- 
Nagumo type reductions. Solitary pulses are constructed 
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in Section III. Section IV studies wave train generation 
at the boundary by periodic firing. In Section V, we de- 
scribe the two main mechanisms for propagation failure, 
related to damage in the myelin sheath and the action 
of chemicals. Section VI contains our conclusions. In 
Appendix A we recall the derivation of the model. Ap- 
pendix B explains the nondimcnsionalization procedure. 
Appendices C, D and E contain additional material on 
pulses. 

II. THE HODGKIN-HUXLEY MODEL FOR 
MYELINATED NERVES 




FIG. 2: Coefficient functions in the HH equations as functions 
of tifc = V, in dimensionless units. 



A. Dimensionless equations 



The dimensionless Hodgkin-Huxley (HH) model for a 
myelinated nerve axon is: 

^ + I{vk,mk,nk,hk) = D{vk+i - 2vk + Wfc-i), (1) 
at 

druk 
dt 

= eA„(wfe)[noo(wfe) ~ Tifc], (3) 

dt 

= eXAhivk) [hoc{vk) - hk] , (4) 



^ Am{vk)[moo{vk) - nik], (2) 



with 



I{v,m,n,h) = gK-n'^iv -Vk) + gNam^h{v - 1) 



(5) 



+gL{v-VL). 

Here, Vk is the ratio of the deviation from rest of the 
membrane potential to a reference potential, Uk is the 
potassium activation, is the sodium activation and hk 
the sodium inactivation. Appendix A recalls the deriva- 
tion of the model. Appendix B details the procedure 



we have followed to nondimensionalize the system. The 
time scale is chosen by looking at the relaxation times 
T„, Tm and Th for Uk, mk and hk- Typically, r„ < t„ 
and r„ ^ Th- Thus, e = ^ is small. 

The rate functions and the stationary states in 
can be fitted to experimental data. Figure [21 plots their 
shape for the motor axon of a frog. Analytic expressions 
for these functions and typical values of the parameters 
for the frog nerve are collected in Appendix B and will be 
used in our numerical tests. The values for gNa, gx , gL 
have orders of magnitude 1, 10^^, 10^^, respectively. The 
coupling parameter D ^ 10^^ and e ~ 10^^ with A ^ 1. 
This means that we have two separate time scales and 
that discreteness effects are relevant. 
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FIG. 3: Temporal action potential Vk{t) (thick solid line) gen- 
erated by a single stimulation at fc = 0. The temporal profiles 
of rrife (dotted), nt (dashed-dotted), hk (dashed) have been 
superimposed, all of them at a fixed node k. The thin solid 
line shows the profile of mt — moo{vk)- All variables are ex- 
pressed in dimensionless units. 

B. Numerical solutions 

System CQl-Q displays excitable behavior when it has 
a unique constant stationary state (v* , m* , n* ,h*), which 
is stable. Figures |21 and 01 show solitary pulses and wave 
train solutions generated by solving (^-Q) numerically 
for the parameter values in Appendix B. After a short 
transient, the system relaxes to a traveling wave: Vk{t) — 
v{k — ct), mk{t) = m(fc — ct), nk{t) ~ n{k — ct) and 
hk{t) = h{k — ct). All nodes undergo the same evolution 
with a time delay i. 

As shown in Figures |31 and 01 pulses and wave trains 
are composed of "sharp" interfaces and "smooth" re- 
gions. Let us describe the temporal profiles of a pulse 
as we move from left to right starting from the equilib- 
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FIG. 4; (a) Same as in Fig. |3] but generated by a periodic 
stimulation at fc = 0. (b) Corresponding spatial profiles at 
a fixed time: Vk (circles), rrifc (asterisks), (squares), ht 
(triangles). Zoom in one spatial period. All variables are 
expressed in dimensionless units. 

rium region (1) in Figure 13 The leading edge of each 
pulse is a sharp interface, marked as (2) in Fig. |3| At 
this leading interface, and /i^ remain almost constant 
whereas Vk and undergo abrupt changes. The evo- 
lution of Vk and rrik is described by a reduced bistable 
system and the leading edge is a wave front solution of 
this system. A smooth region follows, where nk,hk vary 
slowly and Vk and rrifc are quasi-static: — mao{vk) 
and /(ufc, moo(wfe); rife, TOfc) = 0. It is marked as (3) in 
Fig. 121 At the end of this region the approximation 
nik = rriooivk) breaks down, see the thin solid line in Fig- 
ure O A trailing interface develops, where changes 
abruptly whereas nk,hk remain almost constant. This 
new region is marked as (4) in Fig. Notice that the 
fast variable Vk happens to be near its equilibrium value 



V* and decreases smoothly towards it. Diffusion is negli- 
gible and the evolution of Vk and at the trailing edge 
is governed by a system of ordinary differential equations. 
The bistable structure is lost and no trailing wave front 
is formed. 

The variables seem to be split in two blocks: slow 
(evolving in the dimensionless time scale T = et) and 
fast (evolving in the dimensionless time scale t). This 
suggest the possibility of finding an asymptotic descrip- 
tion of pulses by exploiting time scale separation. The 
difficulty of dealing with a discrete space variable is over- 
come by noticing that the traveling wave profiles are 
smooth functions of a continuous variable and solve a 
set of differential-difference equations (see Appendix C). 
Let us first check whether the number of variables can be 
reduced. 



C. Failure of simple FitzHugh-Nagumo type 
reductions 

It is quite tempting to look for a simple asymptotic de- 
scription of pulses in terms of one fast and one slow vari- 

pn 

able, as in the FitzHugh-Nagumo model |3il23l- These 
reduced models usually assume that mk relaxes instanta- 
neously to its equilibrium state = moc{vk) and that 
the sum of the two slow variables is constant during an 
action potential + hk = r 



auie! 
3: 



dvk 
~dt 



+ gKnt{vk-VK)+gNamoo{vk)^{r - nk)ivk-l) 

-HiLivk-VL) ^ D{vk+i -2vk + Vk-i), (6) 



= e A„ (ufe ) [rioo (w/c ) - "-fe] ■ (J) 

In view of the numerical results described in Section fllBI 
these assumptions are inaccurate. The solid line in Fig. 
I^lshows that the difference ~-moo{vk) is not negligible 
at the leading wave front. Setting ruk = moa{vk) distorts 
the speed, as shown in Figure On the other hand, 
nk + hk is not constant during at action potential. We 
may select r = n*+h* to fit the leading edge but a slightly 
different value of r is required for the trailing part of the 
pulse. Keeping Uk + hk ^ r during the action potential 
alters the width of the pulse, as shown in Figure 

In conclusion, the FHN type reduction ©-iQ) may be 
useful to gain insight on pulse motion and propagation 
failure in excitable media, but it is quantitatively inaccu- 
rate for the HH system with realistic parameter values. 

III. ASYMPTOTIC CONSTRUCTION OF 
PULSES 

Accurate descriptions of pulse waves in the HH model 
have to deal with the full set of equations. In this Sec- 
tion, we find an approximation of their temporal profiles 
by matched asymptotic expansions as e ^ 0. This con- 
struction yields predictions for the speed and width of 
the pulses, as well as a characterization of the parameter 
ranges in which propagation fails. We explain the pro- 
cedure for the parameter values indicated in (|B3p . The 
structure of the pulse may change slightly for other pa- 
rameter values. Other possible structures are discussed 
in Appendix IeI 
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FIG. 5: Shape of the source f{v): (a) gNa = 1.49, (b) gNa = 
0.745, (c) gNa ~ 0.497, when (n, h) = {n* , h*) (solid hne) and 
(n, h) — (n'^' , ft'^' ) (dashed hne) . AU variables are expressed 
in dimensionless units. 

In a traveling pulse we distinguish five regions (illus- 
trated in Figure EJ. In each of them, a reduced descrip- 
tion holds. The whole profile is reconstructed by match- 
ing the approximated solutions found in each region at 
zcroth order (see 2^ for a description of this technique). 
Our reference time scale is the slow time scale T = et. 
The technical details of the matching are given in Ap- 
pendix ^ 

The resulting temporal pulse profiles with k fixed have 
the following structure: 

1. Front of the pulse. This is region (1) in Figure |31 
Here Vk ~ v*, mk m* , ^ n* and hk ~ h*, 
(v* ,m* ,n* , h*) being the stable equilibrium state 
of the system. 

2. Leading edge of the pulse, located at T = Tq and 



marked as (2) in Figure O The slow variables, Uk 
and hk, remain essentially constant but the fast 
variables, Vk and mk, undergo rapid changes in the 
time scale t — Xzzla. £ {—oo, oo). To leading order, 



^ + I{vk,mk,nk,hk) = D{vk+i - 2vk +Vk-i), 



dt 



Am(wfc) [moo(wfe) - mk] , 



(8) 



dn 



dt dt 



dvk 



Thus, Uk and nk remain almost constant. Matching 
with the previous region, Uk = n* and hk = h* 
(see Appendix . The evolution of Vk and mk is 
described by the 'fast reduced system': 

I{vk,mk,nk,hk) = D{vk+i - 2vk + Ufc-i), (9) 



dmk 
dt 



^m{vk)[moo{vk) - mk], (10) 



with Uk — n* and hk ~ h* . This system dis- 
plays bistable behavior. Let us denote hy v — 
h), i = 1,2,3 the three solutions of 



f{v;n, h) = -I{v,moo(v),n,h) = 



(11) 



in a neighborhood of {n* , h*). The system (j^- (|10|) 
has a unique stable traveling wave front solution 
Vk{t) = v{k — ct), mk — m{k — ct) joining the 
two stable constant states i'^^'' {n* , h*) , maoii^''^^) 
and i^(^)(n*, h*), m^oii^^^'^), which propagates with 
a definite speed c = c+(n*,/i*). This wave front 
solution is the leading edge of the pulse. 

3. Top of the pulse. This is region (3) in Figure |31 
Here, the slow variables evolve in the time scale 
T and the fast variables relax instantaneously to 
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their equilibrium values: rrik = Tnoo{vk), in which 
Vk solves f{vk',nk,hk) — 0. The evolution of the 
slow variables is governed by the 'slow reduced sys- 
tem': 



equations: 



(12) 



^ = A„ ivk)[noo{vk)-nk], 

^^XAh{vk)[hUvk)-hk], 
for Tq < T < Ti. Matching with the previous stage 
we get Vk = v^^"^ {rik, hk), rik ~ n* and hk = h* 
a.t T — Tq (see Appendix . For the parameter 
values indicated in (jB3p . the third branch of roots 
Vk = v^'^Hnk,hk) of f{vk;nk,hk) = disappears 
at {nk,hk) — (nl-'^l , /il-'^l), colliding with the second 
branch: ^ ^.(2) („[i] ^ /^[i]) = At 

the corresponding time T = Ti, region (3) ends. 
The time Ti is characterized by nfc(Ti) — rS^'^ and 
/ifc(ri) |24|. 

4. Trailing edge oj the pulse, located at Ti and marked 
as (4) in Figure |2| In this region, Vk is no 
longer at equilibrium since the two largest roots 
of f{v;nk,hk) — are lost. The fast variables 
evolve in the time scale t = 21;^ £ (—00,00), 
whereas the slow variables remain essentially con- 
stant. Matching with the previous stage, = n^^\ 
hk = Due to the particular shape of / (see 

the dashed hne in Figure (a)), iy'^^^n^^\ h^^^) is 
close enough to i>'-^'> {n\^\ for Vk — Vk-i to be 
small. Thus, we may neglect the discrete differ- 
ences D{vk+i — 2wfc -|- Wfc-i) and the fast variables 
are governed by a system of ordinary differential 



!^ = -I{vk,mk,nm,hm 



di 



= h.m{vk) [mca{vk) " ruk] 



(13) 



This system has one stable equilibrium point: 
„[2] ^ i,(i)(n[il,/i[il), m[2l = moo(w[2l). The fast 
variables evolve from their initial values Vk = w'^' 
and mfc = m^^l = niooiv^^'^) to the equilibrium 
point. 

5. Pulse tail. This is region (5) in Figure |21 In the 
pulse tail, the fast variables relax instantaneously 
to their equilibrium values: m/j = moo{vk) and 
Vk = v'^^\nk,hk)- The slow variables solve \12\ 
with Vk = v^^"^ (uk, hk) for Ti < T < 00 and evolve 
smoothly from {n^^\h^^^) to their equilibrium val- 
ues {h*,n*) as T ^ 00. 

Uniform approximations to the temporal profiles ob- 
tained by gluing together the approximated solutions in 
each region are given in Appendix IdI Figure compares 
the asymptotic reconstruction to the actual profiles. The 
agreement improves as e decreases. 

Our construction characterizes the speed of the trav- 
eling pulse: it is the speed of the wave front solution of 
©-(Eni with Uk ~ n* and hk ~ h* . The speed of these 
fronts can be predicted by a depinning analysis for small 
speeds (as in ||5|) or in terms of continuous waves for large 
speeds (as in |3j). 

Let us now compute the width of the pulse. In the 



8 



S 0.8 

^ 0.6 
x: 

0.4 
0.2 

Ot 
-0.2. 





■ -A 

























500 



1000 1500 
t 



2000 



2500 



FIG. 6: Asymptotic reconstruction of the temporal profiles 
of (thick solid line), rrik (thick dotted line), (thick 

dashed-dotted line), hk (thick dashed line), versus the actual 
temporal profiles (thin lines). All variables are expressed in 
dimensionless units. 



peak, (rife, lies in the integral curve of: 

~ A„ (n Ji)) [n^ (n Ji)) -n\ 



(14) 



selected by h(n*) = h* and h{n^^^) = n^'^^. Going back to 
the time scale i = — in (|12|l . the peak width is: 



= ^ p; — ^ ^(15) 

' eA„(!^(3)(s,/i(s))) [n„„(!y(3)(s,/i(s)))-sj 



where the integral is calculated along the solution h{n) 
of H14() . Similarly, the duration of tail is: 

ds 



with 



[11 eA„(j.W(s,Ms)))K(i^«(s,Ms))) 



-(16) 



A„(i'(i)(«,/i)) [«^(i/(i)(n,/i))-«] ^ ' 

The integral 72 diverges due to a singularity at n* . How- 
ever, we can use it to predict how long does it take for the 
tail to get close enough to n*, h* replacing n* by n* — rj. 



The spatial length of the peak is found using the trav- 
eling wave structure: Vk{t) — v{k — ct). Thus, Ti is the 
time elapsed from the instant at which the leading front 
reaches a point k to the instant when the end of the peak 
crosses the same point k. The number of nodes Li in the 
pulse peak is approximately the integer part of: 



Li ~c+(n*,/i*)ri. 



(17) 



Our asymptotic construction is consistent when Li > 1. 
This yields a restriction on the size of e for the existence 
of pulses: 

ds 



e < c+ 



-(18) 



/„. K^{vi^){s,h{s)))[n^{ui^){s,h{s)))-s\ 

A similar argument can be applied in the infinite tail 
to quantify the number of nodes at which , hk depart 
noticeably from equilibrium. 

The values predicted by our asymptotic theory for the 
parameters in Appendix B are c = cj^ ~ 0.069, n^^l = 
0.777, /il^l = 0.099, Li = 12, in good agreement with the 
numerical measurements. The leading and trailing edges 
contain about 4 more nodes, that we have neglected. 

Figure [7| compares the predicted and numerically mea- 
sured widths and speeds as the parameters D and g^a 
change. Notice that the values (n*, h*) and {rS^\ /i'^') are 
independent of D. Thus, the dependence of the width on 
D comes through the speed. This explains the similar 
curves observed in Figure [7| (a) and (b). The triangles in 
Figure[7|(a) represent the number of nodes in the leading 
edge, that must be added to compute the total length 
of a peak. Changes in giqa affect the curves v'^'^\n,h)^ 
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FIG. 7: Spatial width of the peak (circles) compared to the 
predicted width (squares) as a function of: (a) D and (c) gjva. 
In (a), triangles represent the number of additional nodes 
forming the leading edge; Speed of the pulse (circles) com- 
pared to the predicted speed (squares) as a function of: (b) 
D and (d) gjva- In (b), triangles represent the speed of the 
FHN reduction. All variables are expressed in dimensionless 
units. 

i = 1,2,3, and the speed of the fronts in the reduced fast 
system. As Figures [T] (c) and (d) show, the quantitative 
agreement between predicted and numerically measured 
widths and speeds is reasonable, as long as we are not 
close to the critical thresholds Z?c and for propagation 
failure. We will discuss this point further in Section Ivl 

For the parameters indicated in HIj3|l , FHN type reduc- 
tions perform poorly. The widths predicted as D changes 
are particularly bad. We find 160 or 250 points in the 
peak, when 18 or 34 are expected. The predictions for 
the speed are better, see the triangles in Figure [T] (b). 
Similar comments apply when gjssa is modified. 



IV. ASYMPTOTIC CONSTRUCTION OF WAVE 
TRAINS 

When a nerve fiber is periodically excited, we expect 
propagation of signals in form of wave trains. Wave 
trains resemble a sequence of identical pulses periodi- 

n 

cally spaced |25j. The asymptotic description of each 
of these pulses is similar to the construction of solitary 
pulses given in Section IIIII However, there are several 
differences. First, the front of the pulses is another pulse 
and not a region at equilibrium. Second, the leading edge 
is a wave front solution of the fast reduced system ©- 
(Unj) with Tife = N, hk = H, {N,H) ^ {n*,h*). This 
wave front solution selects the speed c — {N, H) of 
the wave train. It has to match the tail of the previous 
pulse, described by the slow reduced system ((T^ with 
Vk = (^'^^■'(nfe, hk). Thus, we find a family of wave trains 
for couples (iV, H) lying on an integral curve of: 



dh _ AA^(t/(i)(7i,fe))[/ioo(;^(i)(7i,fe)) -fe] 
d^i ~ A„(z/(i)(n,/i))[noo(i^(i)(n,/^))-n] 



(19) 



The curve is selected by observing that each pulse in the 
wave train has one driving leading edge, as in SectionlTTIl 
This means that /i(n[il) = /il^l, where {n^^\h^^'^) are the 
values for which zy(2)(nW, /jW) = ly^^') {n^^\ h^^^) . 

The spatial period is approximated by adding 
the lengths of the smooth regions: L{N,H) = 
c+{N,H)T{N,n) where T{N,H) = Ti{N,H) + 
T2 {N, H) is the time spent in the excited and recovery 
branches. The temporal lengths 71,72 are defined as in 
SectionHni with n* ,h* replaced by N, H. Now, T2{N, H) 
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of slow and fast variables produces pulses of vanishing 
width. Other scenarios for failure might arise when more 
than two time scales are present. We focus here on mech- 
anisms for failure associated with changes in the param- 
eters due to illness or drugs. 
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FIG. 8: Action potential near propagation failure: (a) Spatial 
profiles of Vk (circles), mjc (asterisks), (squares) and 
(triangles) superimposed at t fixed, (b) Temporal profiles of 
Vk{t) (solid), VTikit) (dotted), nk{t) (dotted-dashed) and hkit) 
(dashed) superimposed for k fixed. All variables are expressed 
in dimensionless units. 

is finite. 

The values predicted by our asymptotic theory for Fig- 
ure H are c = c+ = 0.052, nlil = 0.755, /il^l = 0.09, 
L ^ S7 ~ 38, in good agreement with the numerical mea- 
surements. Numerically, we obtain i ~ 43 for the spatial 
period. The difference is approximately the number of 
points in the edges, that we have neglected. 

V. PROPAGATION FAILURE 

Pulses may fail to propagate at least by two reasons: 
weak coupling and inadequate time scale separation be- 
tween the different variables. Weak couplings cause prop- 
agation failure by pinning the leading edge of the pulse. 
A small time scale separation between the two blocks 



A. Propagation failure due to damage in the 
myelin sheath 



The loss of myelin alters the value of the coupling co- 
efhcient D. Then, the leading wave front can only prop- 
agate if the D is large enough to avoid the pinning in the 
fast reduced system (|^- (|10|) with rife — n* and hk — h* . 
The critical coupling Dc depends on the shape of the 
nonlinear sources, which is controlled by the parameters 
9Na,gK,gL, VNaiVx and Vl- The general rule is that the 
area enclosed by f{v) between its second and the 
third zeroes has to be large enough (depending on D) 
for the leading edge to propagate. The proximity of fail- 
ure is detected by the fact that the wave profiles develop 
'steps'. FigurelHIb) illustrates generation of 'steps' in the 
time profile of the pulse for Vk near propagation failure 
D = 0.0072 ^ Dc- As for bistable equations, the de- 
pinning transition is associated with a bifurcation in the 
system 
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FIG. 9: Impact of decreasing qnu on pulse propagation: (a) 
QNa ~ 0.745, (b) gMa ~ 0.497. The spatial action potentials 
have been superimposed at different times. All variables are 
expressed in dimensionless units. 

B. Propagation failure due to the action of 
chemicals 



uses 
3 



Many drugs block the propagation of nerve impulses 
by reducing sodium and potassium conductivities 
Figure |SJb) illustrates the impact of decreasing gj^a on 
the shape of /. The enclosed area decreases and so does 
the propagation speed, up to a critical value at which 
the width of the pulses (given by formula I^17\i ) vanishes. 
Only decremental pulses are observed, see Figure El 



merical solutions of the discrete HH system. We have 
discussed two mechanisms for propagation failure. The 
first one is related to damage in the myelin sheath and is 
reminiscent of the depinning transitions for wave fronts 
in bistable systems as the coupling decreases. The sec- 
ond mechanism reflects the impact of chemicals on the 
time scale separation in the model. Our analysis applies 
to isolated nerve fibers. However, real motor nerves of 
vertebrates comprise several hundred of interacting fibers 
[2^. It would be interesting to extend our asymptotic 
predictions to bundles of fibers. 

Our asymptotic construction may be useful to under- 
stand systems with a mathematical structure similar to 
HH: models for propagation of impulses through cardiac 
tissue , models of charge transport in semiconductor 
superlattices 2^ or the more complex Frankenhauser- 
Huxley model for myelinated nerves 



piex 
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FIG. 10: (a) Structure of a myelinated nerve fiber (not to 
scale), (b) Equivalent circuit for a myelinated nerve fiber. 

APPENDIX A: THE DISCRETE 
HODGKIN-HUXLEY MODEL FOR 
MYELINATED NERVES 

Myelinated nerve fibers are covered almost everywhere 
by an insulating coat of myelin. Only at small active sites 
(nodes of Ranvier) can the membrane function in the nor- 
mal way. Figure lTUT a') illustrates the structure of a myeli- 
nated fiber. The length of the myelin sheath is typically 1 
to 2 mm (close to lOOd where d is the fiber diameter) and 
the width of the nodes is about l/im. The nodes have 
conduction properties similar to the unmyelinated nerve 
membrane, while the myelin has a much higher resis- 
tance and lower capacitance than the axonal membrane 
|ic| . Myelinated nerve fibers can be described by a linear 
diffusion equation which is periodically loaded by the ac- 
tive nodes [ij, l29l l30| . This picture can be simplified by 



lumping the intcrnode capacitance of the myelin together 

n 

with the nodal capacitance 8]. The myelin is considered 
to be a perfect insulator. This leads to the equivalent 



circuit in Figure ITUT b). C and R represent lumped resis- 
tance and capacitance. Vk, Ik and Iion{k) represent the 
membrane potential, intcrnodal current and ionic current 
at the fc-th node. Applying Kirchoff's laws to the circuit 
yields: 

14-1 ~Vk^ RIk, Ik - Ik+i = + honik) (Al) 

at 



Adopting at each node the Hodgkin-Huxley expression 
for the ion current | , we obtain the discrete Hodgkin- 
Huxley model for a myelinated axon: 



C^+I,on{Vk,Mk,Nk,Hk) 
D{Vk+i-2Vk + Vk-i), 



^ = XMAM{Vk)iM^.iVk) - Mfe), 



^ = XNAN{Vk)iN^iVk) - Nk), 



^ = XHAH{Vk){Ho^{Vk) - Hk), 



(A2) 



(A3) 



where the index k designs the fc-th node of the fiber. 
Here, Vk is the deviation from rest of the membrane po- 
tential, Nk is the potassium activation, Alk is the sodium 
activation and Hk the sodium inactivation. The ion cur- 
rent is given by: 



/,o„(F, Af, iV, H) = g^^ApH{V - VNaM) 



(A4) 



+9LiV - VL,R) + gKN\V ~ Vk,B.)- 
The fraction of open K'^ channels is computed as iV|. 
The fraction of open iVa+ channels is approximated by 
M^Hk- The parameters have the following interpreta- 
tion, gi^a ^^'^ 9k a.re the maximum conductance val- 
ues for Na'^ and pathways, respectively, gj^ is a 
constant leakage conductance. The corresponding equi- 
librium potentials are Vnu, Vk and Vl, respectively. 
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Then, V Na.,R = VNa - Vr, Vk.r = Vk ~ Vb, and 
Vl^r — Vl ~ Vri where Vr is the resting poten- 
tial. C is the membrane capacitance. The coefficient 
D — ^^^ ^^^ ^ = where L is the length of the myelin 
sheath between nodes and ri,re the resistances per unit 
length of intracellular and extra-cellular media. 

This model is adequate for the long axons of peripheral 
myelinated nerves. More biological detail can be included 
by adding an equation for the membrane potential V{x,t) 



across the myelin sheath in the internodes 



3: 



c— - 



1 d'^V V 



— , X e {xk,xk+i),t > {A5) 



dT Ti + re d^x r 

Vixk,t) = Vk{t), V{xk+i,t) = Vk+i{t) (A6) 

coupled with (|A3|I and: 

+ I^oniVk, Mk,Nk, Hk) = hit), (A7) 



hit) = 



1 ,dV , . , dV , 



This model produces a good quantitative approximation 
of the conduction velocity for toad axons 18]. Numerical 
simulations of the sensitivity to different parameters (di- 
ameter, nodal area...) produce results in agreement with 

3 3 



experiments 1191, |20||. The discrete model (IA2ll - (IA3ll is 

recovered by assuming that the axial currents along the 
myelin sheath ^{x,t) are constant in each internode. 



Then, ^{x,t) = ^>'+^W-^''W in [xk,Xk+i] with L 



Xk+i-Xk- As aresult, = L(r.+re) (^fc+i-Vfc+T4-i). 
This approximation is reasonable in view of the numeri- 



cal results in 



APPENDIX B: DIMENSIONLESS EQUATIONS 

For our numerical tests we have selected the param- 
eters and coefficient functions of a frog. For the motor 

n 

nerve of a frog, the data in |33| can be fitted by the fol- 
lowing rate and stationary state functions: 



AM(V-)^0.03[ ^,p^^,;toA\^^-i +4 exp(-|:)], 
M.o(T/)=[l + 4exp(^)22||££^]-i^ 

AH(y)= [0.07cxp(^) + ,^^^,Jo.iv)+i ] ' 



'- 0.07(cxp(3-0.iy) + lj " 

AA^(^) = 0-79[ e4!=o°:?^)-i +0-^25 exp(-^)], 
Noo{V)=[l + 0.125 exp(^^^^^a^(3=^lZKil-i 



(Bl) 



80 ' (0.1-o.oiy) J 
Typical values of the parameters 3^ are given below: 



D 


9Na 


9k 


9l 


l/28(A'm)-i 


0.57 ^mho 


0.104 nmho 


0.025 fimho 


Vr 


VNa 


Vk 


Vl 


—75 mV 


47 mV 


— 75 mV 


-75 mV 


C 




Xh 


Xn 


2.6- A.TpF 


127 (tos)-i 


1.76 (ms)~^ 


2 (tos)-i 



We nondimensionalize the model by choosing as new vari- 
ables Vk = -77^'° , t = TXm, ruk = Mfe, Uk = Nk, 

V Na,R 

hk — Hk- The dimensionless equations are: 



^ + 9Knl{vk - Vk) + 9Namlhk{vk - 1)+ 



9L{vk - Vl) = D{vk+i - 2vk + Vk-i), 
^ = A„(ufe)[moo(wfc) - mk], 



(B2) 



18| (see Figure 2 therein). 



dt 



dhk 

dt 



A„A„(ufc)[noo(wfc) - nk], 
Xh^hivk) [/loo(wfe) - hk] . 



Set G = CXm- Then, the dimensionless parameters are 
given by: 
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APPENDIX D: RECONSTRUCTION OF THE 
TEMPORAL PROFILE 



gNa 


gK 


gL 


D 


Vk 


Vl 


\l 






9k 


9l 


D 


Vk.r 


Vl,r 


A/v 




G 


G 


G 


G 


V Na.R 


Vnu.R 


A A/ 


Am 



The new rate functions and stationary states are ob- 
tained from replacing V by vVNa,R- In dimension- 
less units, the parameters for a frog nerve become: 



D 


9Na 


gK 


gL 


Vk 


Vl 


A,. 


An 


0.093 


1.49 


0.27 


0.065 








0.014 


0.015 



(B3) 



In our asymptotic analysis, we choose A„ = e as small 
parameter and write Xh — eA, A = ~ 1. 



APPENDIX C: EQUATIONS FOR THE WAVE 
PROFILES 

The wave profiles and speeds solve an eigenvalue prob- 
lem for a system of differential-difference equations: 

-cv^{z) = D{v{z + 1) - 2v{z) + v{z - 1)) 

-I{v{z),m(z),n{z),h{z)), 
-cm^{z) = A,n{v{z))[moo{v{z)) - m{z)], (CI) 
-cn^(z) = eAn{v{z))[nooiv{z)) ~ n{z)], 
-ch^{z) = e\Kh{v{z))[hoc{v{z)) -h{z)\. 

In a solitary pulse, the profiles tend to the equilibrium 
states as z — > ±oo. In a wave train, the profiles are 
periodic: v{z) — v{z + L), m{z) — m{z + L), n{z) = 
n{z + L) and h{z) — h(z + L), L being the spatial period. 



We describe below the matching conditions and the 
uniform reconstruction of the pulse profiles in the dif- 
ferent regions. The superscripts (I),...,(V) refer to the 
reduced descriptions corresponding to regions (1),...,(5). 
We have: 

• The matching conditions for the reduced descrip- 
tions of the front of the pulse and the leading edge 
at r = To are: 

V* - vi''\^) « 1, m* - )(^) « 1, 

if e <C To — T <C 1 . The uniform approximations of 
the profiles for T < Tq are: 



V 



unif 



m 



unif 



JII)(T~-To' 



unif 
k 



hl'^'f = h*. 



• The matching conditions for the reduced descrip- 



tions of the peak of the pulse and the leading edge 
at T = To are: 

™^(.f)(nr)(T),/.r)(T))) « 1, 

n^"'\T) - n* « 1, h^"'\T) ~h* 

if e <C T — To <C 1 . The uniform approximations of 
the profiles in regions (II)-(III) are: 

,un.f ^ ,(-)(2:^) + ,(3)(4"^)(T),/.f ^)(T)) 



n 



unif 



Mil) 



unif 
k 
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• The matching conditions for the reduced descrip- 
tions of the peak of the pulse and the traihng edge 
at T = Ti are: 

.f)(4-^)(T),4-^)(T))-.r(-^)«l, 

m^(.f)(4-^)(T),4-^)(r))) « 1, 

ni"\T) - n* < 1, 4''^(r) -h*^l, 

if e <C Ti — T <C 1 . The uniform approximations of 
the profiles in regions (II)-(III)-(IV) are: 

-vf\n\h*)-vf\n\'-\h\'-^), 

+m^{.f\nr\T)Mr\T))) 

n™-^ = n['"\T), hl"'^ = h['"\T). 

• The matching conditions for the reduced descrip- 
tions of the tail of the pulse and the trailing edge 
at T = Ti are: 

-^^^("^(T),C'(T))-.r(2:^)«l, 

n'^^\T)^n* ^l, 4^)(r)-/i*<l, 

if e <C T — Ti <C 1 . The uniform approximations of 
the profiles for T >Ti are: 

- ™r ^ ( ^) + (.j^) inZ' (T) , h^P (T))) 



APPENDIX E: PULSES FORMED BY TWO 
WAVEFRONTS 



For the choice of parameters indicated in l|B3(l . the 
third region of the pulse ends when two branches of roots 
of the cubic source collapse. For other parameters, the 
situation might be different. It might happen that along 
the integral curve H14|l we find a couple {'n}-^\ h}-^^) such 
that there are wave front solutions of the reduced fast 
system traveling with speed c^{ny^\h^'^^) = c+{n* ,h*). 
Then, the top of the pulse ends at this point and the 
fourth region is now a traveling wave front. The pulse is 
formed by two rigidly moving traveling wave fronts. 

Whether a traveling wave front is formed in the back 
of the pulse or not, can be guessed from the shape of the 
nonlinear source /. Figure [^a) depicts f{v;n,h) with 
the parameters values (|B3|) when {n,h) = {n*,h*). It 
is strongly asymmetric. The magnitude of the speed of 
the leading wave front is intuitively related to size of the 
area enclosed by /(u; n*, h*) between its second and third 
zeroes. A*. If we vary {n,h) along the curve 114() . the 
speed of a back front is related to the area A enclosed by 
f{v; n, h) between its first and the second zeroes. We find 
that such areas are always smaller than A* and the cubic 
structure is finally lost. A trailing wave front moving 
at the same speed as the leading wave front cannot be 
formed. Back wave fronts can only be observed for more 
symmetrical sources, when varying {n, h) along the curve 
()14|l we can make A equal to A* . 

We describe here how to modify the asymptotic con- 



struction in Sections IIIII and IIVI to account for pulses 
formed by two rigidly moving wave fronts. In the asymp- 
totic description of these pulses we distinguish again five 
regions. The first three are similar: 

• The front of the pulse is described by Vk ''^ v*, 
ruk ^ m* , Uk ^ n* and ^ h* . 

• The leading edge of the pulse is a wave front solu- 
tion of the fast reduced system (|^- H1U|I with Uk = 
n* and hk = h* joining {u^^") {n* ,h*)^mao{v*)) — 
(w*,m*) and {v^^\n* ,h*),moo{v'^^\n* ,h*)). This 
front propagates with a definite speed c = 
c+{n\h*). 

• In the transition between interfaces, Vk = 
i^^^\nk,hk), ruk = TOoo(wfc) and nk,hk solve the 
slow reduced system (|12ll . evolving from {n* , h*) to 

The fourth region is different: 

• The trailing edge is the wave front so- 
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lution for the fast reduced system ®- 

(jlOII with Uk = nl-'^l and hk — h^^\ join- 
ing (:y(3)(„[l]^;,[l])^^^(j,(3)(^[l]^;j[l])) and 

(i/(i)(n[iU''l),"^oo(j^(''("''U''l)), respec- 
tively, n^^l and /i^^l are now selected in 
such a way that this front travels with speed 
c = c_(n[il,/i[il) =c+(n*,/i*). 

The fifth is similar: 

• In the pulse tail, Vk = iy'^'^\nk, hk), rUk = m^{vk) 
and nk,hk solve the slow reduced system H12|l . 
evolving from (n^^l, /i^^l) to {n*,h*) as t — > oo. 

The temporal and spatial width of the peak can be 
computed as in Section IlIII Now, the number of points 
in the leading and trailing wave fronts are neglected. 

For wave trains, the construction in Section IIVI has 
to be modified as follows: {N,H) and {n^^\h^^^) must 
lie in the same integral curve of H19(l and satisfy c — 
c_ (nil], /ill]) = c+{N,H). 
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